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ABSTRACT 
The  inventory  control  problem  is  examined.   A  FORTRAN  program  is 
developed,  which  evaluates  the  performance  of  an  inventory  against 
simulated  demands  derived  from  generators  according  to  Poisson,  nega- 
tive binomial,  and  normal  probability  laws.   The  program  is  tested 
using  data  obtained  from  an  actual  inventory  control  point  and  four 
sets  of  simple  decision  rules.   The  results  indicate  that  simulation 
is  a  good  way  to  test  decision  rules  and  to  demonstrate  the  results  of 
arbitrary  restraints  such  as  budgetary  restrictions. 
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1.    Introduction. 

Not  too  long  ago  "inventory"  and  "wealth"  were  synonyms.   The  wealth 
of  a  medieval  merchant  was  his  inventory,  the  sum  of  all  his  goods  and 
possessions.   Money  was  one  form  of  inventory,  though  a  particularly  use- 
ful one,  since  it  could  be  exchanged  readily  for  other  forms.   Interest 
was  both  illegal  and  immoral. 

In  some  parts  of  the  world,  even  today,  a  marriageable  daughter 
displays  a  sizeable  portion  of  her  father's  wealth  by  wearing  the  inven- 
tory of  clothing  and  jewelry  she  will  take  with  her  as  a  dowry.   The 
jewelry  is  mostly  in  the  form  of  coins  strung  together  to  form  necklaces 
or  bracelets,  and  their  number  and  size  are  a  point  of  family  pride. 

In  the  southwestern  part  of  the  United  States  a  rancher  proudly 
shows  his  cattle  to  guests.   It  is  considered  poor  taste,  however,  to 
enquire  as  to  their  number,  for  this  form  of  inventory  is  a  good  index 
of  his  financial  status. 

Most  American  businessmen,  whether  private  or  government,  have  atti- 
tudes toward  inventory  quite  different  from  those  illustrated  above. 
Business  inventory  is  thought  to  be  at  best  a  necessary  evil  and  at 
worst  "the  graveyard  of  American  business"  Ij-M. 

Inventory  takes  away  the  freedom  to  use  money  to  meet  obligations 
or  to  engage  in  more  profitable  investment.   Inventory  (more  precisely, 
speculation  in  inventory)  is  blamed  for  excesses  in  business  cycles.   In 
government,  inventory  represents  treasury  borrowings  and,  for  this  rea- 
son, is  considered  to  be  subject  to  interest  charges. 

The  purpose  of  inventory  is  to  relieve  supply  from  constant  depend- 
ence on  manufacturing  or  other  commercial  sources.   If  the  customer  is 


willing  to  wait  until  the  item  requested  can  be  obtained  from  commercial 
stocks  or  from  production,  there  is  no  point  in  maintaining  an  inventory, 
except,  of  course,  that  buying  in  quantity  reduces  costs  attributable 
to  purchase  order  processing.   But  that  is  another  story  and  beyond  the 
scope  of  this  paper.  The  purpose  of  inventory  in  the  Navy  is  "to  par- 
tially uncouple  demand  from  commercial  supply."   The  expense  of  inven- 
tory is  the  cost  of  service  to  the  customer.  Given  that  management  of 
the  inventory  is  optimal,  the  question  of  inventory  size  boils  down  to  a 
definition  of  the  degree  of  service  desired. 

In  the  past  the  Navy  has  formulated  budgets  to  maintain  inventories 
expressed  in  terms  of  "months  of  supply".   A  month  of  supply  is  the 
amount  of  merchandise  expected  to  be  withdrawn  from  the  inventory  per 

month,  based  on  historical  experience  and  expressed  in  dollars.  Our 

2 
inventories  are  composed  largely  of  insurance  items.   Expressed  in 

months  of  supply,  they  usually  look  outrageously  large  to  the  civilian 
leadership  in  the  Department  of  Defense  and  in  Congress,  whose  experi- 
ence has  been  largely  in  private  business. 

In  private  business,  the  justification  for  inventory  is  profit.   To 
be  sure,  profit  in  a  competitive  market  depends,  at  least  partly,  upon 
customer  good  will.   To  the  extent  that  good  will  demands,  a  commercial 
resale  activity  maintains  inventory,  even  though  that  inventory  may  be 
uneconomical.   Within  this  restraint,  it  works  strenuously  to  weed  out 

I  am  indebted  to  Cdr.  S.  W.  Blandin,  S.C. ,  U.S.N*.  ,  for  this  ex- 
pression as  well  as  for  many  of  the  ideas  expressed  in  this  paper. 

''An  insurance  item  is  a  slow-moving  or  non-moving  item  which  is 
retained  in  inventory  for  service  or  essentiality  considerations. 


slow-moving  items.  The  penalty  for  low  or  nonexistent  inventory  is 

occasional  loss  of  sales.   In  the  Navy,  the  penalty  may  be  much  higher. 

The  shortage  penalty  has  been  defined  as  the  cost  of  a  com- 
pletely inoperable  ship  per  unit  of  time.   It  is  obvious,  however, 
that  not  all  shortages  will  render  the  ship  completely  inoperable. 
Some  shortages  will  only  partially  degrade  the  effectiveness  of 
the  ship.   Partial  degradation  is  accounted  for  with  the  concept 
of  essentiality;  essentiality  being  a  relative  measure  of  the 
seriousness  of  a  part  shortage. 

Essentiality  is  determined  from  two  basic  elements :  (1)  mission 
effect,  and  (2)  compensability.   "Mission  effect"  refers  to  the 
operational  capability  of  the  system  when  the  part  is  missing,  and 
"compensability"  refers  to  the  ability  to  make  up  for  the  shortage 
by  repair,  local  manufacture,  substitution,  or  cannibalization. 
Compensability  is  primarily  a  matter  of  time;  that  is,  it  deter- 
mines how  long  the  mission  effect  loss  will  be  felt.   Judgement 
of  qualified  personnel  is  used  to  classify  items  with  numerical 
values  for  mission  effect  and  compensability.   Essentiality  is  the 
product  of  these  two  numerical  values. 

The  shortage  cost  associated  with  a  part  shortage  is  the  product 
of  essentiality  and  the  shortage  penalty  [llj  . 

Shortage  penalties  for  representative  ships  are  as  follows  jjll]  : 


Shortage  Penalty 


Ship  Type 

Per  Year 

Per  Month 

Per  Day 

Guided  Missile 
Heavy  Cruiser 

$11,566,000 

$963,000 

$31,600 

Attack  Aircraft 
Carrier  (Forrestal) 

24,712,000 

2,060,000 

67,600 

Radar  Picket 
Destroyer 

2,244,000 

186,500 

6,130 

Radar  Picket 
Nuclear  Submarine 

8,500,000 

708,000 

23,300 

Ocean  Minesweeper 

690,000 

57,500 

1,890 

The  above  figures  are  based  on  depreciated  construction  and  conversion 
costs  and  total  annual  maintenance  and  operating  costs. 

Although  the  concept  of  essentiality  has  been  understood  for  years, 
and  although  it  has  always  been  applied,  on  a  more  or  less  intuitive 


basis,  ve  have  not  yet  been  able  to  identify  a  significant  number  of  the 
items  in  the  Navy  inventory  with  meaningful  essentiality  codes.   We  are, 
therefore,  unable  to  describe  the  penalty  cost  of  low  inventory  in  a 
budget  request,   Essentiality  is  one  approacn  to  more  meaningful  budget- 
ing. 

Apart  from  essentiality,  Navy  inventory  managers  need  to  justify 
budgets  in  terms  of  effectiveness,  the  ability  to  satisfy  anticipated 
demands.   We  need  to  explain  to  Department  of  Defense  officials  and  to 
Congress  our  need  for  inventory  dollars  in  terras  of  performance,  not 
months  of  supply  or  inventory  to  sales  ratios.   We  need,  most  of  all, 
some  way  to  demonstrate  to  ourselves  and  to  others  the  expected  results 
of  budget  decisions. 

In  the  Navy  as  in  any  large  organization,  we  often  develop  complex 
procedures  and  install  them  without  adequate  testing.   Such  procedures 
frequently  fail  because  they  are  based  on  an  erroneous  or  partial 
apprehension  of  reality.  Or,  on  the  other  hand,  the  procedure  may  be 
correct  enough  and  still  fail  because  of  faulty  execution.   Sometimes 
the  wherewithal  to  make  a  procedure  work  is  missing,  whether  it  be  funds, 
personnel,  time,  computer  capacity,  or  some  other  essential  ingredient. 

As  an  illustration  of  both  kinds  of  failure,  I  offer  the  attempt  to 
identify  recurring  and  non-recurring  demand.   A  demand  which  is  non- 
recurring at  a  retail  issue  point  may  well  be  recurring  when  considered 
system-wide.   For  example,  suppose  a  particular  equipment  is  overhauled 
at  a  shipyard.  As  far  as  the  shipyard  is  concerned,  there  is  a  one- shot 
requirement  for  repair  parts.   But  the  same  equipment  will  be  overhauled 
next  quarter,  only  on  a  different  ship  at  a  different  shipyard.   Clearly 


the  repair  parts  requirements  are  recurring,  yet  it  is  impossible  to 
recognize  them  as  such  at  the  retail  level,   But  even  more  significantly, 
the  definitions  of  "recurring"  and  "non-recurring"  are  so  unclear  and 
subject  to  various  interpretation  that  any  conclusion  drawn  from  such 
categorization  is  wishful  thinking,  particularly  when  we  remember  that 
it  is  frequently  personnel  at  the  GS-3  level  who  make  the  distinction. 

We  are  all  familiar  with  the  failure  of  some  elaborate  system  which, 
having  been  installed  with  great  fanfare,  later  fails  decisively.   We 
are  fortunate  when  the  failure  is  recognized,  and  the  procedure  is 
scrapped.   Normally  what  happens  is  that  the  procedure  is  kept  because 
we  are  "committed"  or  because  we  cannot  bring  ourselves  to  admit  the 
mistake. 

What  is  required  is  some  way  to  evaluate  decisions  before  we  are 
totally  committed.   It  is  the  purpose  of  this  thesis  to  show  how  a 
simple  simulator  can  help  both  in  budgeting  and  in  the  evaluation  of 
inventory  control  stockage  rules. 


2.    Inventory  Policy. 

The  inventory  control  problem  consists  of  two  major  questions: 
when  to  reorder  and  how  much.   In  the  ideal  situation  the  demand  to  be 
experienced  during  lead  tirae  is  known ,  as  is  the  lead  time.   By  lead 
time,  we  raean  the  total  elapsed  tirae  between  recognition  of  the  need  to 
reorder  and  actual  receipt  of  the  merchandise.   The  reorder  point  is  set 
at  precisely  the  known  demand  corresponding  to  the  known  lead  tirae.   The 
order  quantity  is  set  by  a  formula  which  balances  ordering  costs  against 
holding  costs  and  calculates  an  "economic  order  quantity'',   AC  the  ex- 
piration of  lead  tirae  an  optimal  quantity  resupply  is  delivered  just  as 
the  last  unit  of  the  old  inventory  is  shipped   Customer  service  has 
been  perfect,  and  total  expense  has  been  minimized. 

Life  in  the  real  world  of  inventory  management  is  not  so  predict- 
able as  in  the  Utopia  described  above.   Demand  fluctuates  widely  and 
unpredictably;  lead  time  varies;  inventory  restrictions  arc  arbitrarily 
imposed  through  budget  restraints.   Management  reacts  to  these  realities 
by  striving  for  the  best  balance  between  service  and  inventory  invest- 
ment.  Policies  are  translated  into  decision  rules  which  control  the 
answers  to  the  questions,  when  to  reorder  and  how  much. 

The  U.  S.  Navy  Ships  Parts  Control  Center  attains  budget  objectives 
by  adjusting  two  arbitrary  constants  and  by  restraining  economic  order 
quantity  and  risk.   Let  us  examine  SPCC's  system  more  closely  as  an 


^Much  of  the  advanced  thinking  on  Navy  inventory  management  has 
been  done  at  the  Ships  Parts  Control  Center,  Mechanicsburg.  Pennsylvania. 
The  material  in  this  section  is  drawn  largely  from  the  reports  written 
by  the  Special  Assistants  for  Advanced  Logistics  Research  and  Develop- 
ment, of  S.P.C.C. 


illustration  of  scientific  inventory  management  in  the  real  world. 

The  formula  for  economic  order  quantity  used  at  SPCC  is  as  follows: 


EOQ  =  "\  4Q  x  C  x  2  ■  J 

\|    I  x  P      ^l 

where : 

Q  =   quarterly  demand 

C  r  order  cost 

I  s  holding  cost  rate 

P  r  unit  price 

"}     z     a  constant  used  for  attaining  management  objectives 

(which,  being  interpreted,  means  living  within  the 

budget). 

SPCC  imposes  a  further  control  by  restricting  the  economic  order  quantity 

to  values  between  one  quarter's  demand  and  five  years'  demand. 

At  SPCC  risk  is  defined  as  the  acceptable  probability  of  running 

out  during  lead  time.   Risk  is  calculated  as  follows: 

Pxl  x  EOQ    n 
R  =   4Q  x  s   x  A  a 

where : 

Pxl  =  holding  cost 

4Q  =  annual  demand 

S  =  shortage  cost 

EOQ  z  economic  order  quantity 

yl  s  a  multiplier 

Risk  is  normally  restricted  to  values  between  one  percent  and  50  percent. 
Protection  level,  which  is  defined  as  one  minus  risk,  is  used  in  the 


computation  of  variable  safety  level. 

At  this  stage,  we  must  leave  statistics  and  begin  to  use  the  tech- 
niques of  probability.   No  matter  how  stable  our  usage  rates  night  be  - 
and  usage  rates  in  the  Navy  are  almost  never  stable  -  any  forecast  of 
future  usage  is  sure  to  be  somewhat  in  error.   For  thi  reason,  we  base 
our  reorder  point  on  some  type  of  projection  of  past  demand  £lus  a 
safety  level  to  allow  for  variability.   The  amount  of  this  variable 
safety  level  depends  on  the  protection  we  desire  for  each  line  item.   If 
we  eliminate  safety  level  and  set  reorder  point  at  exactly  the  forecast 
of  demand  expected  during  lead  time  we  run  a  fifty  percent  chance  of 
running  out  if  demands  are  normally  distributed  about  some  mean  or 
average  value. 

The  average  demand  is  forecasted  using  exponential  smoothing,  a 
technique  originated  by  Robert  G.  Brown  QJ .   The  formula  is  as  follows: 
Q  =  aD  +  (1-a)  Q 
where: 

Q  s  new  forecast  of  quarterly  demand 

a  =  smoothing  weight  or  constant 

D  =  demand  experienced  during  past  quarter 

Q  =  forecast  for  past  quarter 
The  beauty  of  exponential  smoothing  is  that  it  retains  the  effect  of 
past  experience  without  our  having  to  store  large  amounts  of  raw  data, 
Another  desirable  quality  is  that  later  data  are  weighted  more  heavily 
than  older  data,  the  value  of  the  weighting  being  a  function  of  the 
smoothing  constant.   A  further  convenience  lies  in  our  ability  to  vary 
the  smoothing  constant  to  obtain  either  gradual  or  rapid  response  to 
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change,  as  desired.   For  a  new  item,  for  example,  the  constant  can  be 
set  high,  say  at  .3  to  .5.   As  usage  data  are  acquired,  the  constant  can 
be  lowered,  to  .1  perhaps,  or  even  Lower,  so  as  to  give  less  relative 
effect  to  last  quarter's  demand  and  more  to  older  history. 

We  are  interested  not  only  in  the  average  demand  experienced  in  the 
past  but  also  in  the  variations  from  this  average.   Specifically  we  are 
interested  in  determining  whether  or  not  past  demands  fit  any  known 
probability  distribution.   If  they  do,  we  can  use  this  distribution  to 
predict  the  probability  of  our  being  able  to  satisfy  forecasted  demand 
with  a  given  stock. 

By  careful  analysis,  SPCC  las  found  that  its  items  fit  three  distri- 
butions fairly  well.   If  the  quarterly  demand  is  .5  or  less,  SPCC  assumes 
a  Poisson  distribution.   If  the  demand  is  greater  than  .5  but  less  than 
or  equal  to  25,  the  negative  binomial  distribution  is  assumed.   If  the 
demand  is  greater  than  25,  it  is  assumed  to  be  normally  distributed. 
These  distributions  are  considered  further  in  Appendix  I. 

To  see  how  the  reorder  point  includes  both  forecasted  demand  and 
variable  safety  level  let  us  look  at  a  normally  distributed  item. 

The  protection  level  is  computed  by  recursively  figuring  the  proba- 
bility that  demand  during  lead  time  will  be  less  than  or  equal  to  an 
increasing  number  of  units  of  stock.  When  this  probability  builds  up  to 
the  previously  calculated  protection  level  the  process  is  stopped  and 
the  number  of  units  corresponding  to  that  probability  is  the  reorder 
point  less  program  stock  and  obligations. 

We  see,  therefore,  that  the  inventory  manager  implements  a  budget 
restraint  by  scaling  down  the  protection  level  through  the  adjustment 


of  the  parameters  which  govern  its  calculation.   The  effects  of  such  a 
restraint  are  not  felt  immediately,  for  lead  times  for  our  technical 
items  are  in  the  order  of  quarters,  or  even  years.   The  effect  of  a 
lower  budget  for  stock  replenishment  will  begin  to  manifest  itself 
toward  the  end  of  the  smallest  lead  time  and  will  become  increasingly 
evident  thereafter. 


Forecasted  Demand 


Safety 
Level 


Units  of  Stock — ^ 


Protection  Level 


The  forecasted  demand  protects  the  lower  half  of  the  demand  distribution. 
The  addition  of  variable  safety  level  takes  us  up  to  the  protection  level 
previously  calculated.  The  area  to  the  right  of  the  protection  level  is 
the  risk. 
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Reorder  point,  then,  is  calculated  as  follows: 
RP  s  LQ  +   VSL  +  P  +•  0 
where : 

L  r   lead  time  in  quarters  plus  one  standard  deviation  to 

protect  against  lead  time  variations 
Q  s  quarterly  demand  forecast 
P  s  program  stock 
0  z     obligations  (back  orders) 
VSL  =  variable  safety  level 
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3.    Inventory  Control  Simulation. 

In  1962,  International  Business  Machines  Corporation  completed  its 
Modular  Inventory  Management  Simulator,  an  elaborate  computer  program 
in  the  FORTRAN  compiler  language.   The  program  is  segmented  so  as  to 
permit  great  flexibility.   It  is  possible,  for  example,  to  use  one  of 
the  forecasting  routines  in  the  program  or  to  insert  the  routine  used 
loci  1  y.  The  program  includes  all  the  functions  normally  associated 
with  inventory  management:  demand  forecasting,  error  measurement,  sales, 
orders,  receipts,  etc. 

The  program  written  for  this  thesis  is  an  attempt,  on  a  much  smaller 
scale,  to  show  how  demand  can  be  simulated  according  to  probability  dis- 
tributions and  how  such  simulated  demand  can  be  used  to  evaluate  inven- 
tory performance  under  various  sets  of  decision  rules.  The  demand  gener- 
ators are  subroutines,  which  can  be  used  in  any  program  where  random 
variable  generators  are  required.  The  program  itself  is  reproduced  in 
Appendix  II. 

The  first  ten  instructions  of  the  main  program  accomplish  the  pre- 
liminary housekeeping  tasks  of  recording  the  number  of  line  items  under 
consideration,  setting  aside  storage  space  for  the  variables,  zeroizing 
the  distribution  counters  and  starting  the  random  and  normal  (0,1) 
number  generators. 

The  second  segment  reads  in  the  raw  data,  computes  the  probability 
distribution  parameters  for  each  item,  identifies  the  distribution 
applicable  to  each  item,  counts  the  number  of  items  applicable  to  each 
distribution  and  calculates  for  each  item  the  average  quantity  per 
requisition.   In  a  normal  situation  these  data  would  be  continuously 
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available  as  one  of  the  requirements  of  normal  business,, 

The  expected  demand  for  each  item  is  calculated  by  averaging  the 
first  three  quarters  and  then  applying  the  exponential  smoothing  formula 
to  this  average  and  the  fourth  quarter  demand.  A  high  smoothing  con- 
stant, 0.3,  is  used  because  the  data  came  from  a  new  inventory  manage- 
ment center  and  its  demands  are  not  believed  to  have  settled  down  yet- 
Similarly,  for  the  first  three  quarters  the  squared  demands  are  averaged. 
The  expected  squared  demand  is  computed  from  this  average  and  the  fourth 
quarter  squared  demand  by  means  of  exponential  smoothing  with  the  same 
smoothing  weight,  0.3.  The  variance  is  then  obtained  by  the  familiar 
formula: 

cr  2  =  e(x2)  -tj<*3 

In  assigning  items  to  probability  distributions  it  is  necessary  to 
test  the  standard  deviation- to-mean  ratio  of  items  having  average  quar- 
terly demands  of  over  25.   If  this  ratio  exceeds  three  the  item  is 
identified  as  negative  binomial  distributed.  This  step  is  required 
because  the  probability  of  having  a  negative  demand  must  be  zero,  or 
practically  zero.   If  the  mean  is  three  times  the  standard  deviation  and 
the  normal  distribution  holds,  the  probability  of  a  negative  demand  is 
negligible  (.0035). 

The  program  now  completes  setting  up  the  problem  by  recording  stock 
quantities  in  accordance  with  decision  rules  which  may  be  simple  or  very 
elaborate.  These  rules  are  in  a  subroutine,  which  can  be  changed  at 
will.  The  inventory  is  priced,  extended,  and  totaled. 

Finally,  the  program  secures  simulated  demands  from  the  demand  gen- 
erators, extends  and  totals  issues,  and  calculates  effectiveness,  the 
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ratio  of  total  quantity  issued  on  demand  to  total  quantity  demanded. 

FORTRAN  is  the  computer  language  used  for  this  thesis.  More  spe- 
cifically, the  program  was  written  for  the  FORTRAN  60  compiler  for  the 
Control  Data  Corporation  1604  computer. 

FORTRAN  permits  the  use  of  symbolic  machine  coding  interspersed 
with  FORTRAN  statements.  Symbolic  coding  is  used  in  the  subroutines 
which  generate  pseudo  random  numbers  and  pseudo  random  normal  (0,1)  ran- 
dom variables;  and  its  use,  though  highly  desirable,  does  reduce  the 
generality  of  the  program,  for  the  symbolic  language  is  that  of  the  1604. 
Other  machines  have  different  symbolic  languages,  and  the  instructions 
would  require  translation  before  the  program  could  be  run  on  another 
computer.  One  possible  solution  would  be  to  substitute  generators 
written  for  another  machine.   Such  generators  are  widely  used  and  should 
be  obtainable. 

A  second  limitation  in  these  two  subroutines  is  that  both  presume  a 
48-bit  word  length. 

The  Control  Data  1604  general  purpose  digital  computer  at  the  Naval 
Postgraduate  School  was  available  for  this  project.  This  computer  has 
a  storage  capacity  of  32,768  48- bit  words. 
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4.   Demand  Generators. 

The  simplest  of  the  three  demand  generators  is  that  for  the  normal 
distribution.  A  subroutine  was  available  at  the  Postgraduate  School 
which  generates  pseudo  random  normal  (0,1)  variables.  The  number  secured 
from  this  subroutine  is  converted  by  multiplication  by  the  standard  devi- 
ation of  the  item  and  then  adding  the  predicted  quarterly  demand.  The 
result  is  rounded  off  to  the  nearest  whole  number. 

For  the  Poisson  and  negative  binomial  items  a  pseudo  random  number 
in  the  range  zero  to  one  is  secured  from  a  generator,  which  was  also 
available.  This  number  is  regarded  as  a  cumulative  probability  mass. 
The  probability  masses  for  the  item  under  consideration  are  calculated 
and  summed  thus:  P(O)  -f-  P(l)  -)-P(2)  -f-. . , .  until  the  total  mass  equals 
or  exceeds  the  random  number.  The  last  and  highest  value  of  the  random 
variable  is  the  simulated  demand.  Computationally,  it  is  simpler  to 
compute  P(0),  subtract  the  random  number,  and  then  continue  to  calculate 
and  add  probabilities  to  this  difference  until  it  becomes  either  zero  or 
positive. 

For  the  Poisson  distribution 

P(0)  =  e"*» 

where : 

P(0)  *  the  probability  that  the  random  variable  is  zero 
e  s  the  natural  base,  2.718 
q  s  the  forecasted  quarterly  demand 
Thereafter,  probabilities  are  calculated  recursively  according  to  the 
formula 

P(D   s   p^i-1)  x  q 
i 
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where ; 

P(I)   r   the  probability  that  the  random  variable  is  equal 

to  I 

P(I-l)  r  the  probability  that  the  random  variable  is  equal 

to  1-1 

q  r   the  forecasted  quarterly  demand 

The  iterations  are  terminated  when  I  a  10,  since  the  probability  that  the 

random  variable  will  exceed  10  is  negligible. 

For  the  negative  binomial  generator  the  parameters  p,q,  and  r  must 

be  computed  as  follows: 

p  •  AM 
V 

q   =   1-p 

r   :   AM  X  p 

q 

where ; 

p  a  the  probability  of  success  for  a  single  Bernoulli  trial 

q  s  the  probability  of  failure  for  a  single  Bernoulli  trial 

r  s  the  number  of  successes  (called  AR  in  the  program) 

AM  «  the  forecasted  quarterly  demand 

V  ■  the  variance  in  quarterly  demand 

The  random  variable  is  the  number  of  failures  before  the  rth  success. 
The  probability  of  zero  is  thus : 

P(0)   s  pr  ,  and  thereafter 

P(I)   »  P<I-1)  x  (r  +  I-l)  x  q 
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In  the  iterative  scheme  of  this  program  the  largest:  value  that  is 
permissible  for  I  is  16,383,  this  figure  being  the  upper  limit  for 
indexing  a  "do  loop".  This  limitation  could  be  overcome  at  the  expense 
of  a  few  more  instructions. 

The  purpose  of  all  the  demand  subroutines  is  to  generate  random 
"noise"  around  the  demand  forecasted  by  some  manipulation  of  historical 
data. 
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5.   Testing  the  Program. 

The  raw  data  used  for  the  test  consisted  of  four  quarters  of 
demand  information  on  277  items  in  Federal  Supply  Class  2815,  Diesel 
Engines  and  Repair  Parts.  The  data  were  secured  from  the  Defense 
Construction  Supply  Center,  Columbus,  Ohio,  and  are  typical  of  demand 
in  repair  parts  classes  in  that  high  variance~to-mean  ratios  are  preva- 
lent. The  breakdown  into  fast,  medium,  and  slow  movers  was  not  typical 
since  94  (34;,)  of  the  items  had  average  demands  of  over  25  per  quarter, 
155  (567)  had  average  quarterly  demands  of  25  or  less  but  greater  than 
or  equal  to  0.5,  while  only  28  (101.)  had  average  quarterly  demands  of 
less  than  0.5.  The  corresponding  (estimated)  percentages  on  ail  cur- 
rent SPCC  items  are,  by  way  of  contrast,  14.8%,  26.17.,  and  59.1%, 
respectively. 

It  is  recognized  that  277  items  and  four  quarters  of  usage  data 
are  not  enough  for  statistical  accuracy.  They  are  enough,  however,  to 
show  whether  or  not  the  program  works. 

When  the  data  were  smoothed  and  checked  for  high  variance- to-mean 
ratios,  only  eight  of  the  277  iteras  were  identified  as  normally  dis- 
tributed, while  245  were  negative  binomial  distributed,  and  24  were 
Poisson  distributed. 
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Four  sets  of  very  simple  stockage  rules  were  tested. 


Forecasted  Quarterly 

Demand 

Unit  Price 

Stock 

First  Rules 

Poisson  Items 

-^  .25 

>-  $100 

2  Ea 

^  .25 

^  $100 

3  Ea 

i£  .25 

^  $100 

1  Ea 

£.  .25 

^  $100 

2  Ea 

Negative  Binomial 

and 

Normal  Items 

All 

-^  $100 

U 

All 

$100^U/P-^$10 

U  +  S 

All 

^  $10 

U-I-2S 

Second  Rules 

Poisson  Items  Same 

as 

Above 

Negative  Binomial  and  Normal  Items 

Ail  All  U 

Third  Rules 

Poisson  Items  Same  as  Above 

Negative  Binomial  and  Normal  Items 

AH  -^  $50  U 

All  ^  $50  u-f-S 

Fourth  Rules 
Poisson  Items 

^  .25 

Otherwise 

Negative  Binomial  and  Normal  Items 
All 
All 
All 
All 

Where  U  is  the  forecasted  quarterly  demand  and  S  is  the  quarterly 
demand  standard  deviation. 
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«-==-  $100 

1  Ea 

0 

^  $100 

.5U 

$100^U/P-^$50 

U 

$50^U/P^$10 

u-fs 

^  $10 

U+2S 

For  each  set  of  stockage  rales,  stock  levels  were  established.   50 
simulation  runs  were  then  made  using  demands  from  the  demand  generators. 
The  average  results  are  shown  below.  They  show  how  effectiveness,  dollar 
sales,  and  sales- to- inventory  ratios  are  affected  by  changes  in  stockage 
decision  rules.  The  evaluation  of  these  and  any  other  results  depends 
on  the  goals  of  management. 

In  the  first  three  sets  of  rules,  the  24  Poisson  items  were  held  at 
relatively  high  inventory  levels.  For  the  remaining  items,  various  com- 
binations of  unit  price  criteria  were  used  to  determine  the  size  of  the 
variable  safety  level.  The  second  set  of  rules  is  particularly  inter- 
esting  in  this  respect  because  it  shows  the  result  to  be  expected  if  no 
safety  level  is  provided. 

The  fourth  set  of  rules  drastically  reduced  the  inventory  levels 
for  Poisson  items  and  applied  more  elaborate  criteria  to  the  remainder. 
The  result  was  a  dollar  inventory  of  only  76  percent  of  the  first  one, 
which  yielded  comparable  effectiveness  and  a  better  sales -  to- inventory 
ratio.  This  result  was  possible  because  the  seldom  demanded  Poisson 
items  tended  to  be  high  priced. 


Inventory 

Sales 

Sal 

es/Inven 

tory 

Effectiveness 

First  Rules 

$276,615.98 

$143,501.35 

.518 

,952 

Second  Rules 

$182,009.56 

$118,389.19 

.651 

.680 

Third  Rules 

$239,530.38 

$136,101.93 

.569 

.878 

Fourth  Rules 

$211,419.04 

$118,329.89 

.561 

.947 

6.  Conclusions. 

In  significant  positions  -  notably  at  the  top  -  the  Department  of 
Defense  is  staffed  with  civilian  personnel  who  understand  mathematical 
models,  gaming,  and  probability  theory.   These  people  rightly  insist  on 
discussing  budget  problems  in  such  terms.   It  behooves  us  in  the  military 
to  learn  to  speak  their  language.   Inventory  control  simulation  with 
probabilistic  demand  generation  is  one  way  to  do  so. 

To  be  useful  to  an  inventory  control  point,  a  simulator  needs  to  be 
more  elaborate  than  the  one  written  for  this  thesis.   It  should  be  con- 
sistent with  the  other  routines  of  the  center,  so  that  it  would  be  pos- 
sible, at  any  time,  to  leave  off  normal  processing  and  project  as  far 
into  the  future  as  is  desired,  using  simulated  demand  and,  if  desired, 
revised  decision  rules.   It  should  take  into  consideration  both  antici- 
pated receipts  (dues)  and  back  orders  (obligations).   It  should  establish 
dues  and  obligations  and  should  record  receipts  and  the  release  of  obli- 
gations. 

The  IBM  Modulator  Inventory  Simulator  does  all  these  things,  To- 
gether with  probabilistic  demand  generators,  it  could  form  the  basis  of 
a  local  program. 

The  Navy  is  always  faced  with  the  task,  not  only  of  understanding 
its  problems,  but  of  communicating  them  to  the  civilian  leadership  and 
to  Congress.  Inventory  control  simulation  can  help  us  do  both.   It  is, 
first  of  all,  a  good  way  to  look  at  our  own  business  and  at  proposed 
changes.   It  is,  further,  a  way  to  present  budgets,  showing  clearly  the 
anticipated  results  of  various  levels  of  funding.   It  deserves  serious 
consideration  as  a  management  tool. 
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APPENDIX  I 

PROBABILITY  DISTRIBUTIONS 

The  normal  probability  law  is  too  well  known  to  require  treatment 

in  this  paper.   It  is  thoroughly  discussed  in  almost  any  probability 

text. 

Less  familiar  is  the  Poisson  distribution,  also  called  the  law  of 

rare  events.   It  is  defined  as: 

The  number  of  occurrences  of  events  of  a  specified  type  in  a 
period  of  time  of  length  1  when  events  of  this  type  are  occurring 
randomly  at  a  mean  rate  y(_    per  unit  time  [ij . 

The  Poisson  probability  mass  function  is 

e    A.     x  =  0,1,... 
x; 

Its  mean  and  variance  are  both  equal  to  A.. 

Even  less  familiar  is  the  negative  binomial  or  Pascal  distribution, 

This  distribution  is  defined  as: 

The  number  of  failures  encountered  in  a  sequence  of  independent 
Bernoulli  trials  (with  probability  p  of  success  at  each  trial) 
before  the  rth  success  [Vj . 

The  negative  binomial  mass  function  is 


( 


r+  X- 1\   r  v  s        ^ 

)  p  qx    x  =  0,1,...  .    r  =  1,2,...  .    0£  p±   1  . 


x 


Its  mean  is  — °     while  the  variance  is  -r~?£  .   The  negative  binomial 
P  p2 

fits  the  demand  distribution  of  many  items  of  supply.   One  of  its  desir- 
able properties  is  that,  while  the  probability  of  the  random  variable 
being  equal  to  zero  may  be  quite  high,  it  is  impossible  for  it  to  be 
negative.   This  is  a  characteristic  which  makes  the  negative  binomial  law 
fit  demand  experience  of  items  having  relatively  high  variance  to  mean 
ratios,  where  the  normal  distribution  would  be  inappropriate . 
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Main  Program 


is   1*1 


Set  N 

i 


For  each  item 
derive 

paramo ters   & 
forecast  usage 


Establish 
inventory 


Extend   & 
total 

■     — 


USETOT, 
DEMSXT,  ANIJM, 
DENOM 


Get 
DEMAND(I) 


T 


DEMAND(I)  x 
PRICE(I) 

— *  DEMEXT 


ANUM«*DEMAND(I) 

>*  ANUM 

DEMAND(I)  x 
PRICE(I) 
— 5>  USEXT(I) 


ANUM-*STOCK(I) 

— **  ANUM 
STOCK(I)  x 
PRICE(I) 

— *-  USEXT(I) 


USETOT  +  USEXT( 

>-  USETOT 

DENOM  -t- 
DEMAND(I) 

5NOM 


I) 


anum/deno> 

iFFECl 

I 1 
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END 


Poisson  Demand  Generator 


Get  R 


e-^PROB 


N  +  1-*N 


*& 


PROB-R 


SUM 


JlJL£_ 


1-*N 


M-^AN 


0— >AN 


probxq/an 


PR  OB 


I 


SUM+PROE 
->3UM 


End 
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Negative  Binomial  Demand  Generator 


1  +  1-*I 


Compute  P, 

Q,  &AR 
Round  off  AR 


Get  R 


»  PROB 


PROB  -  R 
— *  SUM 


Yes 


1-*I 


I-^AN 


0-*AN 


(AR-f-AN-1)  x 
PROB  x  Q  /  AN 
— >  PROB 


SUM  «f  PROB 
— *  PROB 


No 
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Normal  Demand  Generator 


Get  DEV 


3  x  DEV  -p  AM 

*>     AN 


Round   off  AN 


End 
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PROGRAM  DATA1 

N=?77 

DIMENSION    QUSAGEt  300 ) , VAR ( 300) ,REQAVE  (300) , 
1  ISTKNRt  300) ,ATIME( 300 ) , S I GMA ( 300 ) , 
2INDIC(300)  .PRICEt'300)  ,  STOCK (300), 
3DEMAN0(300),SrKEXT(300),USEXT(300) 

COMMON    QUSAGE,     STOCK,     INDIC,     PRICE,    SIGMA 

NUNIF=1695 

IPINOM=0 

NEG3IN=0 

IPOIS=0 

CALL    RANCOM(0,R) 

SETUP    PECULIAR    TD     TEST    DATA 

DO    15    1=1, N 
10READ    2    ISTKNRt I ),PRICE(I ) , I REC2 . I QREC 2 . 
1NREC2,IQNREC2, IRE C3 , I QRECi, NREC3, IQNREC3, 
2IREC4,  IQRECU.NREIU, IQNREC4 , IREC 1 , IQREC1, 
3NREC1, IQNREC1 
2    FOR MAT ( I7,12X,F7. 2,5X,8( 12,14) ) 

NR1     =     IREC1    +    NREC1 

'NR2    =    1RLC2    +    NREC2 

NR^    =    IREC3    +    NREC3 

NRU    =    IRtCU    ♦    NREC4 

QTY1=IQREC1+IQNREC1 

QTY2=IGREC2-UGNREC2 

QTY3=IQREC3+IQNREC3 

QTY4=IQREC4+IGNREC4 

NR    =    NR1     +    NR2    +    NR3    +    NRU 

QTY  =  QTY1+QTY2+QT\T3+QTYU 

QUSAGEt I )=(QTY2+3TY3+QTY4)/3. 

AVESQ=(QIY2»QIY2*QTY3»QTY3+QTY4*QTY4)  /3. 

ALPHA=.3 

QUSAGEt  I  )=ALPHA»3TYH(  1 . -ALPHA  )  *QUS  AG  E  (  I) 

AVESQ=ALPHA»OTYl»QTYl+( 1 . -ALPHA ) *AVES Q 

CALCULATION    OF    V<*RIANCE    AND    STANDARD    DEVIATION 

VAR(  I  )=AVESQ-QUSAGE( I  )*QUSAGE( I ) 

SIGMAd  )=SQRTF(  VAR(  I  )  ) 

IDENTIFICATION  OF  DEMAND  DISTRIBUTION! 

IFtSIGMAt I )-0USA3Et I )/3. )  8,8,4 

4,4,3 


IF  (QUSAGEt I)-25. ) 

INDIC { I )=1 

IBINOM=IBINOMH 

GO  TO  7 

IFtQUSAGEt I)-. 5)  6,6,5 

INDICtI )=2 

NEGBIN=NEG8IN+1 

GO  TO  7 

INDICtI )=3 

IP0IS=IP0IS+1 

ANR  =  NR 

REQAVEt I )=QTY/AN* 

STKSUM=0. 

CALL  BUYKN) 

EXTEND  INVENTORY  AND  TOTAL 

0041     I»1.N 

STKEXTt  I  {-PRICE (I  )* STOCK (I) 

STKSUM=STKSUM+STKEXT{ I ) 
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JfcT  DEMAND  FROM  DEMAND  GENERATORS  AND  EXTEND 
USETOT=0. 
DEMEXT=0. 
ANUM=0. 
DENOM=0. 
0048^  1  =  1  >N 

IF( INCIC( I)-2)  42,43,44 

42  CALL  NORMAL  ( NUN! F , QUSAGE ( I ) , S I GMA ( I ) , AN) 
GO  TO  45 

43  CALL  NEGBINOIQUSAGEU  ),VAR(  I),  AN) 
GO  TO  45 

44  CALL  POISSON(QUS<VGE(  I  >,  AN) 

45  D£MAN0( I ) =  AN 
DEMEXT=DEMCXT+DEMAND{ I)*PRICE( I ) 
CALCULATE  EFFECTIVENESS 
1F(ST0CK(  I  )-DEMAMD( I ) )  49,46,46 

49  ANUM=ANUM+STOCK(I ) 

USEXTU  )=  STOCK (  I)»PRICE(  I  ) 
GO  TO  47 

46  ANUM=ANUM  +  DEMAND(  I ) 
•USEXTU  )=DEMAND(I  )* PRICE!  I) 

47  USETOT=USETOT+USEXT{ I ) 

48  OENOM=DENOM+DEMA^D( I ) 
EFFECT=ANUM/DENOM 

52    PRINT    51    STKSUM,JSETOT, EFFECT 
51    FORMAT( 1X,3(F15.3) ) 
END 


Ol 
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C      POISSCN  CtMAND  GENERATOR 
SUBROUTINE  PO  ISSDN ( Q, AN) 

1  CALL  RANDOM! 1,R) 

2  PROB=2.718*»(-Q) 

3  SUM=PROB-R 
U  IF(SUM)  7,5,5 

5  AN=0. 

6  GO  TO  13 

7  00  12  N=l,10 

8  AN  =  N 

9  PROH=PROB»Q/AN 

10  SUM=SUM+PROB 

11  IF(SUM)12,13,13 

12  CONTINUE 

13  END 
C 

C  NEGATIVE    BINOMIA.     DEMAND    GENERATOR 

SUBROUTINE    NEGBHO    (AM,    V,    AN) 

1  P  =  AM  /  V 

C      P  RESTRICTED  TO  VALUES  BETWEEN  .00000  1  AND  .999999 
;IF(P-. 000001)  202,2.200 

200  IF(P-. 999999)  2,2,26l 

201  P=. 999999 
GO  TO  2 

202\P=. 000001 

2  0  =  1  .-P 

3  AR=AM«P/C+.5 

5  IAR=AR 

6  AR=IAR 

7  IF(AR)  8,8,10 

8  AR=1. 

10  CALL  RANDOM  { 1 , R) 

11  PROB  =  P**AR 

12  SUM  =  PROB  -  R 

13  IF(SUM)  16,1U,1U 
1U  AN=0. 

15  GO  TO  22 

16  DO  21  1=1,16383 

17  AN=I 

18  PROB=( AR+AN-1 . ) A N»PROB»Q 

19  SUM  =  SUM  +  PROB 

20  IF(SUM)  21,22,22 

21  CONTINUE 

22  END 
C 
C      NORMAL  DEMAND  GENERATOR 

SUBROUTINE  NORMAL  ( NUNIF , AM, S, AN) 
1  CALL  RNDEV  (NUNIF.DEV) 
3  AN=S«DEV*AM+.5 

N=AN  < 

v         AN  =  N 

U  END  m 
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C      NORMAL  (C,l)  RANDOM  VARIABLE  GENERATOR 

SUPROUTINE    RNUEVt NUN1F.DEV) 

CON<NZ=7777777777777777B) . 

RSV{ ERA=1  ). 

ENA(O)  STACERA) 

SILKS)  FN  II  (  15) 

8RDV    LDA(NUNIF)  MUK  1220703125) 

SCL( 7600000000000 OOOB)       STA(NUNIF) 

RAO(ERA)  1  J  PI  (8R0V) 

SUFH200000000000DOOOB)       STA(ERA) 

AJP2CL+1)  LAC(ERA) 

ENAKC)  LLS(36) 

SCAK2057B)       LRS(U7) 

SSK(ERA)  SCM(NZ) 

SCM(NZ)  FMU(  172467331727205UB) 

8    STA(DEV)  LNI1 <N). 

RETURN 

END 
C 

C      RANDOM  NUMBER  GEMERATOR 
'SUBROUTINE  RANDOM  (K,R) 
1  IF(K)  5,2,5 
B    2  E  =  UOOO  0000  COOO  0000 
B    3  E2=20C0  COOO  GOOD  0000 
B    U  X  =  1777  7777  7777  7777 

5  SIU6(L+7),ENI6(2) ,ENQ(N) , LDA ( X ) ,LLS ( 1 9), 
SCL(E),ADD(X),SC  (  E  ) ,  IJP6(  L-  1  )  ,  ST  A  (  X)  . 
ARS( 1 1),ADD(E2),FAD(E2),STA(R)tENI6(N), 
ENI(N) 

6  END 


•1 


- 
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C      FIRST  SET  OF  STOCKACE  RULES 
SUBROUTINE  BUY1 [V ) 
10DIMENSICN  GUSAGEl  300)  , STOCK ( 300), 
1  INDIC (30C) , PRICE!  3CC)  ,  SI  CM A (300) 
COMMCN  CLSAGE.ST3CK, INDIC , PR  ICE , S IGMA 
H020  I=1,N 
C      IDENTIFICATION  OF  POISSON  DISTRIBUTED  ITEMS 

IF(  INCICU  )-2)  13,10,2 
C      POISSCN  RULES 

2  IF(GUSAGb(  D-.25)  6,6,3 

3  IF(PRICE(I )-100.)  5,5,U 
U    ST0CK(I)=2. 

GO  TO  20 

5  ST0CK(I)=3. 
GO  TC  20 

6  IFCPW ICEI I I-10C.I  8,8,7 

7  ST0CK(D  =  1. 
GO  TC  20 

8  ST0CK(I)=2. 
GO  TO  20 

C     'NORMAL  AND  BINOMIAL  RULES 
10  IF(  PRICE!  D-100.)  12,12,11 
;  11  ISTCCK  =  CLSAGE(  I  »♦  .5 
GO  TO  19 

12  IF(PRICE(  I  )-l0.  )  1l4,lU,13 

13  ISTCCK  =  GLSAGE(  I  )+SIGMAU  )*.5 
GO  TO  19 

1U  ISTOCK  =  QUSAGE(I )f 2 . »S IGMA ( I)  ♦  .  5 
C      ROUNC  OFF  TO  NEAREST  WHOLE  NUMBER 

19  STCCKU  )  =  ISTOCK 

20  CONTINUE 
25  END 

C 

C      SECONC  SET  OF  ST3CKAGE  RULES 
SUBROUTINE  BUY1 M ) 
10DIMENSICN  QUSAGEI  300) ,STOCK(300), 
1 INDIC (30C), PR ICE(  30C) ,  SIGMA(3CC) 
COMMON  GUSAGE,ST3CK, INDIC , PR  ICE , SIGMA 
0020  1  =  1, N 

C      IDENTIFICATION  OF  POISSON  DISTRIBUTED  ITEMS 
J  IF( INC IC(  I  )-2)  10,10,2 

C      POISSCN  RULES 

2  IF(QUSAGE(I )-.25)  6,6,3 

3  IF(PRICE(  D-100.)  5,5,U 
U  ST0CK(I)=2. 

GO  TO  20 

5  ST0CK(D=3. 
GO  TO  20 

6  IF(PRICE(  D-100.)  8,8,7 

7  ST0CK(I)=1. 
GO  TO  20 

8  ST0CK(D=2. 
GO  TO  20 

C      NORMAL  AND  BINOMIAL  RULES 

10  ISTCCK=CUSAGE(I)+.5 
C      ROUND  OFF  TC  NEAREST  WHOLE  NUMBER 

17  STOCK(I)  -    ISTOCK 

20  CONTINUE 

25  END 


i 
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THIRD  SET 
SURROUTIN 
10CIMENSICN 
1  INPIC(30C 
COMMCN  GU 
0020  1=1, 
IDENTIFIC 
I  F  (  I N  C  I  C  ( 
POISSCN  R 

2  IF(QUSAGE 

3  IF(PRICE( 
U  STOCK!  I  )  = 

GO  TO  20 

5  STOCK(  I  )  = 
GO  TO  20 

6  IFIPfUCEl 

7  STOCK! I )  = 
•  GO  TO  20 

8  STOCK(I)= 
GO  TO  2  0 

4  NORMAL  AN 

10  IF(PRICE( 

11  ISTOCK=CU 
GO  TO  17 
ISTCCK=QL 
ROUND  OFF 
STOCK! I ) 
CONTINUE 
END 

FOURTH  SE 
SUBRCLTIN 
ODIMENSICN 
1 INPIC130C 
COMMON  CL 
C020  1=1, 
IDENTIFIC 
IF( INCIC! 
POISSCN  R 
IF(PR  ICE( 

2  IFIGUSAGE 

3  STOCKt I  )  = 
GO  TO  20 

6  STOCK(I)= 
GO  TO  20 
NORMAL  AN 

10  IF(PRICE( 

11  IST0CK=.5 
GO  TO  19 

12  IF(PRICE( 

13  ISTOCK=QU 
GO  TO  19 

1»4  IF(PRICE( 

15  ISTOCK=0U 
GO  TO  19 

16  ISTOCK=QU 

19  STOCK(I)= 

20  CONTINUE 
25  END 


OF  sto: 

E  BUY1 (N 
QUSAGE! 
),PR.ICEl 
SAGE,ST3 
N 

AT  ION  OF 
D-2)  10 
ULES 
(  I)-. 25) 
D-1G0.) 
2. 

3. 

I  1-100.) 
1. 


KAGE  RULES 

) 

3CC) ,ST0CK(3CC) , 

300) ,SIGMA(300) 

CK,INDIC, PRICE, SIGMA 

POISSON  DISTRIBUTED  ITEMS 
tlC, 2 

6,6,3 
5,5,4 


8,8,7 


2. 

D  BINOMI 
I  )-5C.  ) 
SAGE( I)+ 


12 

17 
20 
25 


1 


AL  RULES 
11.11,12 
SIGMAtI  )+.5 

.5 


SAGE!  I  )♦ 

TO    NEAREST    WHOLE    NUMBER 
=    ISTOOC 


T  OF  STJ 
E  BUYKVJ 
GUSAGE! 
)  , PRICE! 
SAGE, STO 
N 

ATION  OF 
I)-2)  10 
ULES 
D-100.) 
(  I  )-.25> 
1. 

0. 

D  BINCMI 
D-100.) 
•QUSAGE! 

D-5C.  ) 
SAGE! I)* 

D-10.  ) 
SAGE!  DO- 


CKAGE   RULES 

) 

300) ,ST0CK(300), 

300) ,SIGMA(300) 

CK, INDICt PRICE, SIGMA 

POISSON    DISTRIBUTED     ITEMS 
,  10,2 

2,6,6 
6,6,3 


AL    RULES 

12,  12,11 
I  )  +  .5 

lU,1i|f  13 


16,16,15 
SIGMAU  )  +  .5 


SAGE(D*2.»SIGMA(I)  +  ,5 
ISTOCK 
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